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1 Introduction 



Dynamics of self-interacting scalar fields is still a subject of very interesting and 
important research. Because such fields are ubiquitous in physics, from condensed 
matter to cosmology, the results can have a wide impact on our understanding of 
the material world. Nonlinearity of pertinent field equations is the reason for 
the existence of very nontrivial and often unexpected phenomena. In particular, 
much attention has been paid to so called oscillons, extremely long-lived scalar 
field configurations which in many respects behave like solitons, but which finally 
decay into radiation, see e.g. recent papers HI. The nonlinearity can have various 
forms. In the current literature dominate models with polynomial Lagrangians. 
However, in a recent paper [2j we have pointed out a rather interesting class of 
scalar field models with non polynomial field potentials which are V-shaped at 
their minima and which have another kind of nonlinearity. The simplest example 
of such models is the signum-Gordon model. 

The signum-Gordon model (s-G) has the following Lagrangian J3l 

L= l d<p d<p . , (1) 
2 dx^ dx^ 

where ip is a real (l+l)-dimensional scalar field, /i = 0, 1, and x°, x 1 are dimen- 
sionless variables obtained by the appropriate choice of units for the physical time 
and position coordinates, g > is a dimensionless coupling constant. The perti- 
nent field potential has the form 

U((p) = g M- 

Its left and right derivatives at cp = 0, where it has the absolute minimum, do 
not vanish - they are equal to ^fg. Hence, in case of the s-G model the second 
derivative U"(0) is (in a sense) infinite. In majority of field theoretic models the 
field potential U((p) is smooth at its absolute minimum at ip = tp , and A = 
([/"((^q)) -1 / 2 is recognized as the fundamental length scale in the model. In the 
case of massless models A is infinite, while finite A > characterizes massive 
models with the finite mass parameter ml = 1/Aq. In this classification, the s-G 
model corresponds to Ao = and infinite ml - for this reason it can be called a 
supermassive model. 

Several rather interesting features of the models with V-shaped field potentials 
have been pointed out in our previous papers: the existence of static, compact 
solitons (topological compactons) flU , the presence of a scaling symmetry [0, 
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and the variety of exact self-similar solutions (31 0. In these papers we have also 
presented a sound physical motivation for considering such models. It includes 
the description of an infinite chain of harmonically coupled classical pendulums, 
and the dynamics of an elastic string (a vortex) pinned with a constant force to 
a line. Because of such a down-to-earth physics behind the s-G model, it is not 
surprising that the model is perfectly well-behaved from the physical viewpoint. 
In particular, the conserved energy is bounded from below, and the field equation 
is of the standard hyperbolic type with the usual causality features. Nevertheless, 
the dynamics of the scalar field in this model is strongly influenced by the fact that 
the field potential is not smooth at the minimum. Moreover, the model cannot be 
linearized even if the amplitude of the scalar field is arbitrarily small, contrary to 
the models with smooth field potentials. Perhaps it is one of the reasons for which 
the models with V-shaped field potentials were not discussed in literature. 

In the present paper we show that the s-G model possesses solutions with fi- 
nite energy which are periodic in time - the oscillons. These oscillons are rather 
interesting for the following reasons. First, they are given explicitly as exact solu- 
tions of the s-G field equation. Moreover, they do not loose their energy because 
they do not emit any radiation. They have a compact support in the space, i.e., a 
strictly finite size. In these respects they differ from the oscillons of the </> 4 model 
which are known only as approximate numerical solutions and seem to contain a 
radiative component. Because of the scaling invariance of the s-G equation in fact 
there is a continuous family of oscillons with arbitrary positive energies. One can 
also trivially combine such oscillons to form multi-oscillon solutions. 



2 The signum-Gordon equation 

The signum-Gordon equation is the Euler-Lagrange equation corresponding to 
Lagrangian (1). It can be written in the following form 

d 2 tp(x,t) d 2 tp(x,t) , . . 

where t = yfgx , x = i/gx 1 . The sign function has the values ±1 when <p ^ and 
for ip = 0. Because of the — sign(^) term Eq. (1) is nonlinear in a rather special 
way. Mathematical aspects of this equation are discussed in [0. Generally speak- 
ing, physically relevant are solutions which are continuous functions of x and t, 
but only piecewise smooth. They belong to the class of so called weak solutions 
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of partial differential equations. The reader not familiar with the mathematics of 
such solutions should consult the literature, e.g., (6l|71- 

The novel feature of Eq. (2) is that the nonlinear term — sign(t^) remains finite 
for arbitrarily small \cp\ provided that cp ^ 0. This fact has profound influence 
on the dynamics of weak fields. In other models nonlinear terms vanish when 
fields become arbitrarily small, cf. the — ip 3 nonlinear term in the case of the non- 
linear Klein-Gordon equation. Consider, for example, the Gaussian wave packet 
ip(x,t = 0) = exp(— x 2 ) at the initial time t = 0. Let us assume for simplicity 
that d t (p(x,t = 0) = 0. The field acceleration d' 2 (p(x,t = 0) is obtained from 
the field equations. In the case of the — cp 3 nonlinearity the dominating at large 
|x| contribution to the acceleration comes from the (9^(x, t = 0) term and it is 
positive, hence the Gaussian wave packet will spread out. In the case of the s-G 
equation the — sign(y9) term dominates at large \x\ and it is negative, hence the 
wave packet will shrink for certain time. Heuristically, both 'forces' —ip 3 and 
— sign(y?) push the field toward the equilibrium value if = 0, but the — sign(^) 
force is much more effective in this respect. 

The signum term in Eq. (2) remains constant until ip becomes equal to zero. 
Therefore, on each interval on the x axis such that cp has a constant sign on it, one 
can use the well-known formula for the general solution of the one dimensional 
wave equation, suitably modified to incorporate the constant +1 or —1 term in the 
equation: 

<p(x,t) = h(x - t) +g(x + t) + cit 2 , (3) 

where 2ci = ±1(= —sign((p)). The functions h, g and the constant c x are deter- 
mined from the initial and boundary conditions. 

3 The oscillon 

Motivated by the heuristic considerations about the wave packets we shall try to 
construct a solution of the s-G equation (2) which does not spread out at least 
during a certain time interval. We start by specifying simple initial data at t = 0: 

[ if x < 0, 

ip(x, 0) = 0, d t ip(x, 0) = I v(x) if < x < 1, (4) 

[ if x > 1, 

where the initial field velocity v (x) is assumed to be negative. Its exact form will 
be determined later. Moreover, we impose the following boundary conditions 

<p(x = 0, t) = 0, <p{x = 1, t) = 0. (5) 
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The hope is that the field will remain localized in the interval < x < 1 for 
arbitrary times t > 0. 

Because v(x) is negative, for small enough t > the field (p will be negative 
or equal to zero. Therefore, in the region < x < 1 we use formula (3) in which 
ci = 1/2. The functions h and g are determined from conditions (4), (5). For a 
certain reason which will become clear later we denote the solution given below 
by ip~. Simple calculations give the following results 



Lp (x,t) 



where 



(pi(x,t) 
ip 2 (x,t) 
<Ps(x,t) 



(fi(x, t) if < x < t, 
(p 2 (x, t) if t < x < 1 — t, 
(f3(x, t) if 1 — t < X < 1, 



x 1 f t + x 1 . , 

-— +tx + - / dsv(s), 
I 2 Jt-x 



(6) 



1 rt+x 

- / dsv{s), 

I Jx-t 



2 2 j.v-i 

x 2 1 1 r 2 -*-^, . . 

— — — - — xt + x + t + - as vis). 

2 2 2 Jx-t 



(7) 
(8) 
(9) 



The reason for splitting the unit interval < x < 1 into the three parts is 
that the boundary conditions (5) modify the field at the points which are causally 
connected with the points x = and x — 1: such points form the two subintervals 
< x < t, 1 — t < x < 1. These subintervals meet at the time t = 1/2. 
Hence, the solution given above is valid for t in the interval [0, 1/2], and for x in 
the interval [0,1]. 

Note that at the initial time t — only (p 2 is present because the spatial sup- 
ports of </?i, (fz are shrank to the points x — 0, x — 1, respectively. The functions 
(fi, (f 2 match each other at the point x = t, similarly as (p 2 , y?3 at x = 1 — t, so 
that <p~ as well as d t <p~ , d x (p~ are continuous functions of x, t. 

In the next step we extend the solution ip~ to the whole x axis. Of course 
we just would like to put ip(x, t) = for x < and for x > 1. This is possible 
provided that 

d x i Pl (x = 0,t) = 0, d x ip 3 (x = l,t) = 0, (10) 

because otherwise c^y? would have discontinuities at the points x = 0, x — 1. 
Such discontinuities, if present, would have to move with 'the velocity of light' 
±1 because the s-G equation belongs to the class of hyperbolic partial differential 
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equations. Conditions (10) give the following equations 

v(t) = -t, v(l-t) = -t, 

where < t < 1/2. It follows that the initial field velocity v(x) has the following 
form 

r — x if < x < \, 



v ^-{ x -l if \<1<1, 



see also Fig. 1. 




(ID 



Figure 1: The initial profile of d t ip for the oscillon solution. 



Now we can compute the integrals present in formulas (7-9). It is convenient 
to split the time interval [0, 1/2] into two parts. It turns out that when t has values 
from the interval [0, 1/4] then 



(pi(x,t) 
<p 2 (x,t) 

<p 3 (x,t) 



x 

T' 

t 2 



— xt if 

| + t 2 - f - | + | if 

\ + t(x-l) ' if 
{l-xf 



t<x<\-t, 
\-t<x<\ + t, 
\+t<x<l-t, 



(12) 
(13) 
(14) 
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and for t from the interval [1/4, 1/2] 



( *\ J ~T if 0<x< |-t, 

^ = ifU- f -| + | if i-t<x 2 <t, (15) 

= t+ y-2"2 + 8' (6) 

^ rJ) = {-^ if i + t<x<l. (17) 

Note that 9 t v?(x, t) = at the time t = 1/4. 

To summarize, we have found a solution of the s-G equation such that ip = 
outside the interval [0, 1] on the x-axis, and inside this interval it has the form 
given by formulas (6), (12-17), provided that < t < 1/2. Snapshots of the field 
configurations at the times t — 1/8 and t = 1/4 are presented in Figs. 2 and 3. 
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Figure 2: Picture of the oscillon solution at the time t = 1/8. For x from the in- 
tervals [1/8, 3/8] and [5/8, 7/8] the function (p(x, 1/8) is represented by rectilinear 
segments which at the time t — 1/4 shrink to the points 1/4 and 3/4, respectively, 
see the first and third lines of formula (13). 
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A (p (x,t=l/4) 




Figure 3: Picture of the oscillon solution at the time t = 1/4. This configuration 
is 'the turning point' because precisely at this time d t (p changes its sign from 
negative to positive. 

For times t from the interval (1/4, 1/2) the solution has the same shape as pre- 
sented in Fig. 2, but now the two rectilinear segments expand and move until they 
cover the intervals [0, 1/2], [1/2, 1] on the rr-axis. 

In order to find (p for times larger than 1/2 we compute and d t tp~ at t = 1/2 
and we take these functions of x as the new initial data for the s-G equation. The 
calculation is trivial, but one should remember about watching the domains of the 
involved functions because they are given by inequalities which also depend on 
time. The result is quite surprising and crucial for further progress: 

<p-(x,t = 1/2) = 0, dt(p-(x,t = 1/2) = -v(x), (18) 

where v(x) is given by formula (11). These formulas imply that <p(x, t) for t > 
1/2 can be obtained with the help of symmetries of the s-G equation: the time 
translation t — > t — 1/2 and the change of sign of the field <p — > —(p. The solution 
obtained in this way is denoted by (p + . Thus, 

p+( x ,t) = - V -(x,t-^). (19) 

The solution (p + (x, t) holds for t in the interval [1/2, 1]. At the timet = 1 the field 
(p and its time derivative d t (p return to their initial values (4). Thus, the period of 
the oscillon is equal to its spatial size, i.e., to 1. 
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The structure of the oscillon solution during one period is schematically de- 
picted in Fig. 4. The functions y?^, where a = 1A, IB, are just parts of 
formulas (12)-(17) - the new notation is introduced for the clarity of the figure. 
Let us list these functions for the convenience of the reader: 

x 2 _ t 2 t x 1 

^ = -y ^ib = 2 +tX ~2"2 + 8' 

t 2 t 2 x 2 x t 1 

f2A = ~2 " tx ' = 2" + f ( x ~ X )' Pic = f2 + y - g ~~ 2 + 8' 

£ 3/ 3 X. .o 

^3A = 2" + 2 + 2 ~ 8' ^ 3B = ~2^~ ^ ' 
and the functions </?+ (x, t) are obtained from formula (19). Note that the oscillon 
is spatially symmetric: the solution is invariant under the spatial reflection at the 
point x = 1/2, that is with respect to the substitution x — > 1 — x. The point 
x = 1/2 is the center of the oscillon. 




Figure 4: Schematic picture of the oscillon solution. Because of the periodicity in 
time, the figure from the unit square should be repeated for t < and for t > 1. 
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The s-G equation is dilation invariant [3|: if ip(x, t) is a solution of it, than so 
is also 

<pi(x,t) = ZV(-py), 

where / is an arbitrary real, positive constant. Applying this symmetry to the 
oscillon solution presented above we generate a one-parameter family of oscillons 
of the length l/-Jg with the period equal to the length (let us recall that we use the 
dimensionless variables x^). The amplitude of the oscillations of the scalar field is 
equal to Z 2 /16. The total energy of the oscillons is given by the formula obtained 
from Lagrangian (1) 

E = \V9 r dx t(^) 2 + ^) 2 ] +V9 [°° dx M = Z 3 . 

Z J —oo J —oo 

The frequency of the oscillations of the field is ~ I' 1 , hence the energy E vanishes 
in the limit of high frequencies. 

The s-G equation is also invariant with respect to Poincare transformations. In 
particular, boosts provide oscillons moving with an arbitrary constant velocity v 
such that \v\ < 1. 

4 Remarks 

1. Note that the s-G oscillons have a strictly finite size, and that the field ap- 
proaches its vacuum value cp — in the parabolic manner, see the functions 
ViAi ^33- m mese respects they are similar to the static topological compactons 
discussed in [4] in another model with a V-shaped potential. The compactness 
and the parabolic approach to the vacuum field seem to be generic features of the 
models with V-shaped field potentials. 

2. The s-G equation admits explicit solutions which describe arbitrary chains of 
static, non overlapping oscillons of arbitrary sizes. The multi-oscillon solutions 
are obtained trivially by adding appropriately translated in space single oscillon 
solutions. Such oscillons do not interact with each other because they have strictly 
finite sizes. 

3. We have attempted to investigate the stability of the s-G oscillon under small 
perturbations. Our numerical simulations did not show any instability. On the 
analytic side, the problem turns out to be quite difficult. The main reason is that 
there is no linear approximation to the function sign((/?) around ip = 0. This fact 
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seems to render the standard linear stability theory useless. Note however that for 
(p 7^ and the perturbation e small enough 

sign(y? + e) = sign(y?). 

Therefore, in such cases the s-G equation implies the simple massless wave equa- 
tion for the perturbation, 

d^e = 0. 

It follows that the perturbation does not grow: it splits into a left- and right-movers 
which just travel with the velocities ±1 until they reach a point at which the func- 
tion sign((/9 + e) changes its sign. It remains to be investigated what happens at 
such points. Generally, one could expect complicated dynamical processes in- 
cluding an emission of radiation or of small oscillons. 

4. We expect that the s-G oscillons will likely have to be taken into account when 
considering various dynamical processes with the s-G field. The oscillons with 
arbitrarily small energies can probably be emitted in such processes, thus forming 
a kind of 'infrared' cloud consisting of a number of small, compact oscillons. 
It is clear that the dynamics of the s-G field turns out to be very exciting. 
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